SU{5) grand unification on a domain-wall brane from an E'e-invariant action 
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An 5*17(5) grand unification scheme for effective 3 + 1-dimensionaI fields dynamically localised on a 
domain- wall brane is constructed. This is achieved through the confluence of the clash-of-symmetries 
mechanism for symmetry breaking through domain-wall formation, and the Dvali-Shifman gauge- 
boson localisation idea. It requires an Ee gauge-invariant action, yielding a domain-wall solution 
that has Ee broken to differently embedded 50(10) ® U{1) subgroups in the two bulk regions on 
opposite sides of the wall. On the wall itself, the unbroken symmetry is the intersection of the 
two bulk subgroups, and contains SU{5). A 4 + 1-dimensional fermion family in the 27 of E^ 
gives rise to localised left-handed zero-modes in the 5* © 10 ® 1 © 1 representation of SU{5). The 
remaining ten fermion components of the 27 are delocalised exotic states, not appearing in the 
effective 3 -I- 1-dimensional theory on the domain-wall brane. The scheme is compatible with the 
type-2 Randall-Sundrum mechanism for graviton localisation; the single extra dimension is infinite. 

PACS numbers: 



I. INTRODUCTION 

If our universe is a 3 + 1-diniensional brane [l|, S H) S S S B S] existing in a 4 + 1-dimensional spacetimc, then 
the most likely field-theoretic origin for the brane is a scalar-field domain wall (DW) or kink [ij. This generic idea is 
naturally compatible with the type-2 Randall-Sundrum (RS2) mechanism for producing effective 3 -I- 1-d gravity on 
the brane [9^| (see, for example, Refs. 13, 11, 12| for the extension of thin-branc RS2 to a domain- wall brane). The 
challenge is to dynamically localise all the other ingredients necessary for a phenomenologically successful effective 
theory on the brane: gauge bosons, fermions, and Higgs bosons. Various localisation ideas for these disparate classes 
of fields have been recently combined to produce an effective brane theory that is plausibly very similar to the standard 
model 
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The purpose of this paper is twofold. First, we wish to point out a very elegant generic connection between the 
clash-of-symmetries (CoS) mechanism for symmetry breaking throu gh d omain wall formation 0, [l5, 1^ ^'^d 



the Dvali-Shifman (DS) idea for dynamical gauge-boson localisation 22] ■ Second, we use this remarkable confluence 



to construct an explicit scheme that realises an SU(5) gauge-invariant effective theory on the brane. In a sense, it is 
a grand unified extension of the model of Ref. but the way in which the Dvali-Shifman mechanism is realised is 
quite different, and we shall argue that it is in fact conceptually more advanced. Remarkably, this scheme immediately 
produces a reahstic spectrum of locahsed fermion zero-modes [i3] (for a review see [2J1) using the simplest possible 
mechanism. While it is beyond the scope of this paper to write down a complete phenomenologically-acceptable 
domain- wall-brane localised SU{5) theory, we shall conclude with brief remarks about how this could be attempted. 
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The clash-of-symmetries phenomenon 14, 15, 16, 17, 18, 19, 2^ automatically arises when the simple Zi kink is 
extended to a theory with a continuous internal symmetry group G in addition to the discrete symmetry. Taking the 
scalar-field multiplet to be in a non-trivial representation of G, the domain-wall configuration spontaneously breaks 
G in addition to reflecting the disconnected vacuum manifold topology created by the spontaneous breaking of the 
discrete symmetry. Two classes of domain- wall solutions exist: those which respect the same subgroup of G at all 
values of the bulk coordinate y, and those where the unbroken subgroup varies in the bulk. We shall call the first class 
"non-CoS domain walls" , contrasted with the "CoS domain walls" of the second class. Clash-of-symmetries DWs can 
arise when the subgroups respected asymptotically (at y — ±oo) are isomorphic but differently embedded subgroups, 
H and H' . The symmetry group at finite y is typically the intersection H n H', which is of course smaller than both 
H and H'. 

The last observation provides an immediate connection with the Dvali-Shifman proposal for dynamical gauge-boson 
localisation. The DS mechanism, as originally proposed [23, envisaged a domain wall configuration where the full 
group G is restored in the bulk, but broken to H in the wall. The gauge bosons of H propagate on the wall either as 
massless Abelian gauge fields or glueballs formed from non-Abelian gauge fields. In the bulk, all gauge bosons have 
to be incorporated into massive G-glueballs.^ Thus, the massless Abelian gauge fields on the wall have to become 
incorporated into massive glueballs in the bulk, and the energy cost associated with the mass gap then plausibly 



localises them to the wall. This heuristic argument is bolstered by the dual-superconductivity model [26|, |27| for 
the confining bulk: the electric field lines from a source charge in the wall are repelled from the interface with the 
dual-superconducting bulk 3, 29|, just as magnetic field lines are Meissner-repelled from an ordinary superconductor. 
The non-Abelian gauge fields of H are also plausibly localised if the mass of the G-glueballs exceeds the mass of the 
ff-glueballs. 

The fact that the full symmetry G is asymptotically restored is clearly not a necessary condition. In the CoS 
situation, the brane-group H D H' is a subgroup of both H and H' , the unbroken symmetries in the two semi-infinite 
bulk regions. By the DS reasoning, provided H and H' contain confining non-Abelian factors, at least some of the 
gauge bosons oi H D H' will be localised. For a realistic theory, we need the localised gauge bosons to include those 
of the standard model. The model-builder needs to engineer the theory to achieve this effect. While this engineering 
shall be the main concern in the rest of the paper, our first generic point has already been made: the clash-of- 
symmetries automatically gives rise to Dvali-Shifman gauge-boson localisation. This CoS alternative realisation of 
the DS mechanism seems conceptually neater than the original, because it can be achieved using scalars in a single 
irreducible representation of G. The original requires two multiplets: a G-singlet to form a kink, which in turn forces 
a G-multiplet to condense in the core of the wall. 

We shall show that the CoS-DS confluence can naturally produce an SU{5) effective theory on the brane. The basic 
ingredients are G — Eq, with the DW-producing scalar held in the adjoint or 78 representation. The groups H and 
H' will be the differently-embedded maximal subgroups 50(10) i^i U{1)e and 80(10)' U{1)e', respectively. Their 
intersection is SU{5)i»U{l)i»U{l), with SU{5) of course being a subgroup of both 50(10) and 50(10)'. Taking both 
of those as conflning gauge theories in the bulk, the localisation of SU (5) gauge bosons follows from the DS-effect. The 
gauge flelds of C/(l) ® J7(l) are not completely localised. When 4 -f 1-d fermions in the 27 of are Yukawa-coupled to 
the scalar multiplet, we shall show that 3 + 1-d left-chiral zero-modes in the phenomenologically-realistic 5* ® 10© 1 © 1 
representation of SU (5) are localised. The remaining ten fermion components remain 4 -f 1-d, and are thus absent 
from the effective brane-theory. The result that the chiralities of the zero-modes come out to be phenomenologically 
correct is very non-trivial, as we shall explain. 

We review the clash-of-symmetries idea in SecHIl Section [Till describes a warm-up example featuring 50(10) CoS 



^ The Dvali-Shifman mechanism requires a confining 4 -|- 1-dimensional gauge theory to five in the bulk. The issue of confinement in 
4-1-1 dimensions is not completely understood, so the Dvali-Shifman mechanism in that context has the status of being a plausible 
conjecture. There is good lattice gauge theory evidence that pure SU{2) Yang-Mills theory with an ultraviolet cut-off is confining in 
4-1-1 dimensions when the gauge coupling constant exceeds a certain critical value [25j | . It is thus plausible that a variety of 4 -|- 1 
dimensional gauge theories exhibit confinement at sufficiently large coupling strength. 
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FIG. 1: The vacuum manifold of a G® ^2 H model. The two circles schematically depict the disconnected coset spaces G/H 
and {G/H)z- Each point along the G/H circle corresponds to a vacuum |0; g) for some g £ G, with the corresponding situation 
for the (G/ H)z circle. The three broken lines represent possible domain wall configurations, with the endpoints a.t y = ±co on 
various choices of vacua. The dotted line represents a possible non-topological domain wall configuration. The short-dashed 
line represents a non clash-of-symmetries domain wall configuration, while the long-dashed line is a clash-of-symmetries domain 
wall. 

domain walls, and explains why the extension to Ef, constructed in Sec II VI is needed. We conclude in SecfVl 

II. THE CLASH OF SYMMETRIES 

Consider a theory (an action) whose symmetry group is the direct product of a continuous symmetry group G and 
a discrete symmetry Z . It is important that Z is not a subgroup of G. Suppose the global minima of the Higgs 
potential spontaneously break G to subgroup H , and simultaneously break Z to a smaller discrete group. For the 
sake of definiteness, we shall take the Z — Z2 — {1, z s.t. = 1} example in what follows, with the Z2 completely 
broken. 

The vacuum manifold then consists of two disconnected copies of the coset space G/H, with the copies related 
by the spontaneously broken z G Z2. This is an immediate generalisation of the simple Z2 kink situation, where 
the vacuum manifold consists of just two disconnected points related by ^2- Each such point is expanded into the 
non-trivial manifold G/H. We shall call the disconnected pieces G/H and {G/H)^. The Z2 must not be a subgroup 
of G for the two disconnected pieces to exist. 

Let |0) be an element of G/H. By definition, h\0) = |0) for all h £ H. Since the Higgs potential is G-invariant, 
if we apply a transformation g G G/H (that is, a transformation such that g G G but g ^ H) to |0), we obtain a 
degenerate vacuum state |0; .g) = .g|0). By considering all possible g's, these transformations generate the G/H piece 
of the vacuum manifold. Applying the non-identity transformation z Z2 from the discrete group, we obtain the 
discrete image |0)z = z\0) of |0). This image is a point in the other disconnected piece {G/H)z of the full vacuum 
manifold. By applying all possible g G G/H to \0)z, the space {G/H)z is generated. Figure [T] illustrates this situation. 

The degenerate vacua |0) and |0; 5) respect differently embedded but otherwise isomorphic subgroups H and iJg, 
respectively. This is elementary: Let hi,h2 G H such that = G H. Then the conjugates 5/11,2,33"^ respect 
the same multiplication table and hence the set gHg~^ is precisely Hg which is isomorphic to H but a different subset 
of G. If h\0) — |0), then trivially ghg~^g\0) = g\0), which simply says that the conjugated elements preserve the other 
vacuum: {ghg^^)\0; g) = \0;g). Similar statements are true for [G/H)^. 

The boundary conditions at y — ±00 for domain wall configurations are chosen from the vacua. If the chosen vacua 
are either both from G/H, or both from (G/ H)z, then the "domain wall" configurations are not topologically stable: 
they are in the same topological class as any of the spatially-homogeneous vacua |0; g), or respectively |0; g) z, and will 
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dynamically decay to one of these vacua. They may be metastable, depending on the Higgs potential topography,^ 
so while they are of some interest we shall not consider them further in this paper. 

Topologically non-trivial DW configurations have one boundary condition from G/H and the other from (G/H)^. 
Evidently, there is an uncountable infinity of such choices, and thus potentially an uncountable infinity of DW 
solutions, all within the same non-trivial topological class. Figured] illustrates the plethora of choices. This potential 
richness has no analogue for the simple Z2 kink. 

Suppose that the boundary condition at y = — 00 is |0) and a,t y — +00 it is |0)2. Then if /i|0) = |0), it also follows 
that 

h\0), = hz\0) = zh\0) = z|0) = 10);,, (2.1) 

because by assumption the symmetry is G ^ so that gz — zg^z Cz Z, g Cz G and hence hz = zh always. Thus, 
the unbroken symmetry at = —00 is precisely the same set H as &t y — +00. A domain wall configuration that 
interpolates between these vacua is then expected to respect the same subgroup H at all y. This is an example 
of a non-CoS domain wall. Clearly, taking the vacua as any pair \Q]g) and lO;^)^ produces a similar outcome (the 
resulting configuration is nothing more than the g transform of the original one). A non-CoS domain wall is the 
simplest possible generalisation of a Z2 kink for a G-invariant theory. 

However, there is obviously a second, more interesting possibility: if the vacuum is |0) aX y = —00, then the 
vacuum at y = +<x can also be a \Q]g)z for 5 7^ 1. In that case, the subgroups respected asymptotically are 
the differeiitly-embedded but isomorphic groups H and Hg, respectively. This defines a CoS-style domain wall 
li E 0, Q, H, H At finite y, the configuration would be expected to respect the smaller group H H Hg due 
to the fact that the solution has to "reconcile" boundary conditions that have different stability groups that "clash" 

So, there will be an infinite family of non-CoS DWs, trivially related to each other by global transformations 
g G G/H. They all have the same energy density, because the Hamiltonian is invariant under G. The CoS DWs have 
a more complicated spectrum. Consider two configurations, Xiiv) ^-nd X2(2/), with xi interpolating between |0) and 
|0; 5i)z, while X2 interpolates between |0) and |0; 32)2, such that gi ^ 32- Suppose, for the moment, that G is a global 
but not a local symmetry. These two solutions cannot be transformed into each by a global G-transformation, so 
they would be expected to have different energy densities (their configurations trace different paths through the Higgs 
potential topography). As a corollary, the non-CoS solutions should have a different energy density from the CoS 
solutions. All these solutions are in the same topological class, so finite-energy dynamical evolution between them 
is allowed. Hence, the special configurations within that topological class that minimise the energy density will be 
topologically stable. The others should be unstable to decay to the minimum-energy configurations, which play the 
role of "vacua" for the "kink-sector" . This general reasoning cannot tell you which configuration has the minimum 
energy- density: you need to calculate that within a specific model. For example, in the toy model considered in 
Ref. lJ| the sign of a Higgs potential parameter determined whether the non-CoS or a CoS solution was energetically 
favoured. 

Suppose now that G is a gauge symmetry, and again consider the configurations discussed in the previous paragraph, 
together with the specification of vanishing gauge fields Am at the solution-level. The non-CoS solutions remain 
connected through global transformations. Two CoS scalar field configurations, xi and X2, can be written as local 
G-transforms of each other. Suppose that 

X2{y) = U{y)xi{y), (2.2) 
where U{y) is a local-G element. Then the original first solution 

X = Xi{y), Am=0 (2.3) 



^ We dread to use the term "landscape" . 

* From experience, we have found that H n Hg is the usual outcome. The specifics depend on the case considered. Sometimes there is 
enhanced symmetry at y = because some of the scalar multiplet components instantaneously vanish there. This enhancement on a 
set of measure zero has yet to find application, although speculations exist [171. 13011 . 
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is gauge-equivalent to 

X = X2{y). Am = ~-{dMU)U\ (2.4) 
e 

where e is the gauge couphng constant, but it is not gauge-equivalent to 

X^X2{y). ^M = 0, (2.5) 

which is the original second solution. Thus the two solutions x = Xii ^-A/ = s-i^d X = X2j Am = have different 
energy densities, even though the scalar-field portions are related by a local symmetry transformation. Although 
Am = — f (9a/J7)C/^ is a pure-gauge configuration, it contributes to the energy density through the x~ Am interaction 
terms. 

Setting the gauge fields to zero at the solution-level is basically a convenient choice of gauge, one we shall adopt from 
now on. Of course the solutions can be made to look very different by gauge-transforming them, but their physical 
consequences cannot change. This circumstance is no different from the monopole or local-string cases, where again the 
solutions look different in different gauges. Actually, it is no more complicated than the usual homogeneous vacuum 
expectation value (VEV) case. If (x) is a homogeneous VEV, then it can be gauge-transformed to a non-homogeneous 
configuration U {x) (x) but the scalar gradient energy is cancelled by the gauge-field contribution. 

The alert reader may have noticed the following: we have not proven that the minimum-energy DW configuration 
must have a gauge-field sector that is gauge-equivalent to zero. This does indeed appear to be a loose end. We shall 
make the assumption that it is in fact true for the purposes of the rest of this paper. Ultimately, one could uncover its 
hypothetical falsity by a perturbative stability analysis for the DW, but that is well beyond the scope of the present 
investigation. 

Finally, a technical point: The set of Hg contains an uncountable infinity of differently-embedded but isomorphic 
subgroups. However, there is a certain useful sense in which the number of embeddings can be considered finite. Let 
the Cartan subalgebra Qc of G be a certain particular set of generators, corresponding to a particular choice of basis 
for the Lie algebra. If we require that the Cartan subalgebras of two subgroups Hg^ and Hg^ are both subspaces of Qc, 
then the number of distinct embeddings is finite. A familiar example of this concerns the SU{2) subgroups of SU{?>). 
While there are an uncountable infinity of ways of embedding SU{2) in 5?7(3), there are only three embeddings that 
have the SU{2) Cartan subalgebras as subspaces of the given Cartan-subalgebra space of SU{3). These are usually 
called I-spin, U-spin and V-spin. When we say "different embeddings" below, this is what we shall mean.^ 



III. SO{W) WARM-UP EXAMPLE: THE NEED FOR Eg 

We now discuss a. G — 50(10) model that serves both as a warm-up for Eq and explains why the extension to Eq 
is necessary. We shall make use of the O(10)-kink analysis of Ref. [17]. While some recapitulation is necessary for 
completeness, we shall be as brief as possible, and the reader is referred to Ref. [l3| for a detailed discussion. 

Let X be a scalar multiplet in the adjoint representation, the 45, of 5*0(10). The most general quartic Higgs 
potential is 

y = lf^' Tr(x2) + i Ai Tt{xY + \\2 Tr(x*), (3.1) 

where x = laX"^ with the A"'s being matrix representations of the generators in the fundamental of 50(10) while the 
/q,'s are the components of the adjoint multiplet. The matrix x is antisymmetric and transforms as per x ~^ AxA'^ 



^ Note that taking linear combinations of Cartan generators to define different embeddings is in accord with Dynkin's general theory of 
embeddings [3l| . In that formalism, the embedding of an algebra H into a simple or semi-simple algebra Q is fully defined by a mapping 
F from the Cartan subalgebra of "H into the Cartan subalgebra of Q, as per Ha — » F{Ha) = X]a=i PaaGa, where Hey, (a = 1, 2, ... , m) 
and Ga (a = 1, 2, . . . , ra) are the Cartan generators of Ti. and Q, respectively. The matrix (Faa) is the defining matrix of the embedding, 
and two embeddings arc different if their defining matrices are different. 
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where A is an SO{10) fundamental-representation matrix. The parameter /i^ is chosen to be positive since Tr(x^) 
is negative definite. The cubic invariant Tr(x'^) identically vanishes so there is an accidental discrete Z2 symmetry, 
X ^ — X: which shall play the role oi Z. It is not a subgroup of 0(10). 

The global minimisation of such a potential was performed by Li [s^l (see also [s^). Using an 50(10) transformation, 
one may always bring a VEV pattern into the standard form 

X = diag(/i e , /2 e , /s e , /4 e , /5 e), (3.2) 

where the fi are real and 

The five independent fields fi correspond to the five generators in the SO{10) Cartan subalgebra. In this basis, 

2 

v^-fj^yyf + x, lyyn +^x2y:ft- (3.4) 



1=1 \i=l / i=l 



For A2 > 0, the global minima of V are at 



where we define /min = \/ M'^/ (lOAi + A2) and the unbroken subgroup is H — U{5). The values of fi at the minima 
are specified up to a sign that can be chosen independently for each component 

f^ = ±/min. (3.6) 

Different choices for these signs correspond to two features: different embeddings of (7(5) in iS'O(lO) and also a choice 
of which Z2 sector the minimum lies in. 

To explore this further, let us turn to possible domain wall configurations. Suppose that aX y = —00, we choose as 
our boundary condition 

x(-oo) = -/min = -/mindiag(e , e , e , e , e). (3.7) 

This defines a certain U{5) unbroken ed y — —00, and the VEV lies in one of the two disconnected pieces of the vacuum 
manifold. At y = +00, there are three choices that lie in the other piece of the vacuum manifold, disconnected from 
the first by the spontaneously broken Z2: 

{/min = /min diag(e ,£,£,£,£) 

/i'/n' = /mindiag(e,e,e,-e,-e) . (3.8) 
fmin = /mindiag(e, -e, -e, -e, -e) 

(Permutations of the minus signs in the last two of these vacua are just a trivial rearrangement of the representation- 
space and need not be separately considered.) Vacua with an odd number of minus signs on the right-hand side on 
Eq. (|3.8p are continuously connected to x(— 00) by 5*0(10) and shall not be considered as they would give rise to 
non-topological domain walls. 

The three vacua in Eq. p.Sp are invariant under differently-embedded subgroups of 50(10): U{5)i, U{5)2 and 
(7(5)3. The superscripts (5), (3, 2) and (4, 1) denote the numbers of plus and minus signs in the VEVs. But they also 
usefully describe the unbroken symmetry of the domain wall at finite y, respectively 

(7(5), C/(3) 0(7(2) and (7(4)0(7(1), (3.9) 



as we now explain. 
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The ansatz for domain wall configurations that interpolate between the stated boundary conditions is xiu) = 
^(2/)x(~oo) + (;(2/)x(+cxd), where the functions h and g obey self-evident boundary conditions. The first configuration, 
which interpolates between — /min and -l-/min for all components fi{y), breaks 50(10) to U{5)i at all values of y, 
because the relative magnitudes of the components are always the same at a given y. It is a non-CoS domain wall. 
The second configuration has an equal-magnitude 3x3 block (of 2x2 submatrices) , and an equal-magnitude 2x2 
block. The unbroken symmetry is then 

t/(3) (g) U{2) = U{5)i n C/(5)2. (3.10) 

Similarly, the third configuration's block structure leads to J7(4) (g) U{1) — U{5)i n U{5)3. 
The Euler-Lagrange equations 



h + 2A2/f , (3.11) 



with the three types of boundary conditions above may be solved numerically. However, a simple way to prove that 

solutions exist is to consider the Ai = slice through parameter space. The equations can then be solved analytically 
to yield 

/i(2/) = /mintanh(/i2/) Vi (3.12) 

for the first boundary condition choice, 

fiiy) = /mill tanh(/i?/) for i = 1, 2, 3 and fi{y) = /mm for « = 4, 5 (3.13) 
for the second choice, and 

fi{y) = /min tanh(^?;) for i = 1 and /^(y) = /^in for i = 2, 3, 4, 5 (3.14) 



for the third choice. The surface energy densities are in the ratios 5 : 3 : 1 for the first to the third solutions [i7[, so 
Eq. (|3.14p gives the topologically stable configuration. 

From a Dvali-Shifman point of view, this stable configuration has an unbroken S'[/(4) on the brane that is embedded 
in SU (5)i on the y < side of the wall, and SU{5)3 on the y > side. The S'C/(4) gauge bosons are thus localised to the 
wall, if the Dvali-Shifman mechanism is correct, because by assumption both SU{5)i and SU{5)3 are in confinement 
phase in their respective bulk regions. This establishes the connection between clash-of-symmetries and Dvali-Shifman 
by way of an explicit rigorously worked-out solution. It is, however, just a toy model since the symmetry breaking 
pattern is not what is required phenomenologically. 

The second configuration, with 

C/(3) (g) U{2) ^ SU{3) ® SU{2) (g) U{1) ® U{1) (3.15) 

on the brane is closer to what we need for a realistic model. While the analytic solution of Eq. (|3.13p is unstable 
to dynamical evolution to Eq. (|3.14p , it could well be that in another region of Higgs-potential parameter space the 
U{3) (g) U{2) solution is the stable one. This has not been established, but let us suppose it is true. The model then 
still does not quite work, although it comes close. 

It is certainly true that the SU{3)(E)SU{2) factor in Eq. (|3.15p is Dvali-Shifman-localised, because it is a subgroup of 
both SU{5)i (the bulk symmetry for y < 0) and SU{5)2 (the bulk symmetry for y > 0). However, there is a problem 
with the hypercharge gauge boson corresponding to U{1)y- To see this, we need to examine the U{1) generators more 
closely. 

The asymptotic gauge groups are 



C/(5)i = SU{5)i U{l)x, and C/(5)2 = SU{5)2 ® U{l)x, ■ 



(3.16) 
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Denote by Yi the hypercharge generator inside SU{5)i, and Y2 the one inside SU{5)2- The two U{iys in Eq. (|3.15p can 
be taken to be generated either by Yi and Xi , or by I2 and X2 , and each pair can be written as hnear combinations 
of the other pair. Now, either Yi or Y2 can be the physical hypercharge Y. Which one is selected will be an accident 
of spontaneous symmetry breaking. At some scale above the electroweak, the breaking 

U{1)y, ® Uil)x, = U{1)y, ® Uil)x, ^ U{1)y, (3.17) 

with either Y = Yi or Y ^ Y2, will have to take place to produce an effective standard model at low-energies (this 
will require an additional Higgs field) . Suppose F = Yi is spontaneously selected. Then the hypercharge gauge boson 
cannot propagate into the y < bulk, but a component of it will propagate into the y > bulk. The generator 
y = Yi is a linear combination of Y2 and X2, so the hypercharge gauge field is a linear combination of the gauge 
fields of Y2 and X2. But only the Y2 part is unable to propagate into the y > region; the X2 part is immune from 
the Dvali-Shifman effect because it is not confining. After electroweak symmetry breaking, this will imply that both 
the photon and Z*^ will leak into the y > bulk, which is phcnomcnologically ruled out. If Y2 happens to become the 
physical Y, then leakage into y < will occur. 

This structure, with localised gluons and bosons, but semi-delocalised photons and Z°'s, almost works. But 
understanding its pathology also provides the cure: We need to expand the symmetry on the brane to contain a 
full SU{5), with the physical hypercharge identified with one of its generators. Further, this brane-S'?7(5) must be a 
subgroup of confining non-Abelian groups on both sides of the domain wall. These two features automatically arise 
when we upgrade from SO{lQ) to Eg as the symmetry of the action. 



IV. THE Eg DOMAIN- WALL BRANE. 
A. Group theory 

Take a scalar field multiplet x in the 78-dimensional adjoint representation of Eq. In the next subsection, we shall 
analyse the associated Higgs potential and produce domain wall solutions. But for now, we just need to use the fact 
that for a range of parameters the global minima of the potential will induce 

Ee-^ SO {10) (E)U{1), (4.1) 

which is a maximal subgroup. Now consider different embeddings of SO{10) (E) U{1) in Eq.^ We shall show below that 
there is a pair of embeddings, which we shall call simply 5*0(10) U{1)e and S'O(IO)' (g) U{1)e', that is of particular 
interest for model-building.^ The domain wall solution we shall find in the next subsection interpolates between x 
VEVs that break Eq to these different but isomorphic subgroups on opposite sides of the wall. The symmetry on the 
wall is then 

[50(10) (E) U{l)E]n [50(10)' E) U{1)e' ] = SU{5) (E) U{1)e E) U{1)e', (4.2) 

as we shall establish. Since SU{5) C 50(10) n 50(10)', the Dvali-Shifman mechanism localises all the SU{5) gauge 
bosons to the domain wall, including the photon and the Z^. 
Let us look at the group theory in more detail. Under 

Ee ^ SO{10)E)U{1)e ^ [SU{5) E) U{l)x] E) U{1)e, (4.3) 

the fundamental 27-dimensional representation of Eq branches as per 

27 -> 1(4) + 10(-2) + 16(1) 

^ l(M) + [5(2,-2) ^5*(-2.-2)]^ [^(-54) +5*(3,i)^10(-i4)]. (4 4) 



® We mean the finite number of embeddings in the sense of the final para graph of Sec. |lT] 
The second embedding has been used in unified model building |34l. l35l [3^ |37(I . 
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The notation for SO{10) (E) U{1)e representations is D{12E), where D is the dimension of the 50(10) multiplet, and 
the ^(1)^; generator has been normahsed as per 

Tr2r{E') = \. (4.5) 

(We use 12E to make the charges integers for convenience. ) The SUib) ® C/(l)x ® U{1)e notation is I}(4v^^,i2£;) 
with 

T^27{X^) = l, Tr27(i;X)=0. (4.6) 
The second embedding is revealed by considering the Hnear combinations 

X' = ~{X + VTEe), E' = ^{-VTEx + E) (4.7) 

that correspond to 

Ee ^ SO(IO)' «) Uil)E' ^ [SU{5) C/(l)x' ] ® C/(1)e'. (4.8) 
Rewriting the multiplets from £)(4VT5X,i2£;) notation to D^^^^-^, i2E') notation, we see that 



l(0,4) 


= 1(-5,1) 


5(2.-2) 


= 5(2,-2) 


5*(-2 -2) 


- ^(3,1) 


l(-5,l) 


= 1(0,4) 


5*(3,1) 


^ ^(-2,-2) 


lo(-i.i) 





so the 5*'s flip roles as do the singlets. Let us now redundantly denote the multiplets through 



The 10 of 50(10)' is 



-(2-2) ^-(-2.-2) 

whereas the 10 of the original 50(10) was instead formed by 



Similarly, the 16 of 50(10)' consists of 



(4.9) 



(4VT5X',12B')' ^ ^ 



5f2!-2i ® 5/0.-2^ (4-11) 



^(2^-2) ® ^(3,1^' (4-12) 



iKi)®5*(-r'^®io!:S!' (4-13) 



whereas the 16 of the original 50(10) consisted of 

lU^e 5X^12) e^lS- (4.14) 

The 50(10)' singlet is i|q ''/^\ whereas the original 50(10) singlet is ^!>^_^^^y 

Because all higher-dimensional representations of E^ are formed from products of 27's, the feature that some 
SU{5) (El U{1)'^ submultiplets flip when {X,E) {X',E') propagates to all irreducible Eq representations. The 
submultiplets can be packaged in 50(10) (SiU{1)e multiplets, or repackaged into 50(10)' ® U{1)e' multiplets. This 
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establishes, constructively, that the two embeddings exist, and that Eq. (|4.2p is true.^ Note that the additional [/(l)'s 
are there because adjoint configurations cannot rank-reduce. 
Let us repeat this exercise for the adjoint of Eq: 

78 ^ 1(0) + 45(0) + 16(-3) + 16*(3) (4.15) 

-,(0:0) 

^ -^(0,0) 

+ ^ -^(0,0) + -^"(-1,-3) + -^^^(1,3) +^^(0,0) J 

r , (-5,-3) .*(3,-3) 1 n(-l,-3) i 
+ ^-^(5,3) +^(3,-3) +-^"(4,0) J 

+ [l!-5l-3)+ 51:3:3! + 10*^2)] (4.16) 

The flipping of roles is evidently 

1(0,0) (0,0) 
-^(0,0) -^(0,0)' 

in*(-4,0) -,(.*(1,3) 
-■-^(1,3) ^ -'-'^(-4,0)' 

-■-^(-1,-3) ^ -'-"(4,0) 
i(-5,-3) (5,3) „x 

^(5,3) ^ -^(-5,-3)- y^-^') 

The S't/(5) adjoint 24[g'°| is common to both 50(10) embedd ings, as befits its status of being in the intersection of 
the two. 

The two i|q q| muhiplets play important roles. Giving a VEV to the 1(0) in Eq. ((iTTS)) breaks Eq to SO{10)(»U{1)e, 
while a VEV for the second singlet in Eq. (|4.16p breaks Eq to 30(10)' (S)U{1)e' ■ A clash-of-symmetries kink interpolates 
between these two VEVs imposed as boundary conditions. At | j/| < 00, both SU (5) (E) U{1)'^ singlet components of the 
78 have nonzero values, and this is precisely why the configuration breaks Eq to the intersection of the two subgroups. 
To analyse this further, we must consider the dynamics. 

B. Higgs potential and domain- wall solutions 

The adjoint scalar multiplet x shall be represented by 

X = /„X", a = l,...,78 (4.18) 

where X's are matrix representations of the generators for the 27 of Eq, and the /'s are the field components. It 
transforms according to 



X 



UxU'^ (4.19) 



where U is group representation matrix for the 27. We shall only be concerned with two of the seventy-eight fields: 
those associated with [E, E'), equivalently {X, E) or (X', E') depending on what basis we choose for the Lie algebra. 
We thus specialise to 

X^fEE + fxX = fEE + fE'E' (4.20) 

with 

fE^fE + ^, (4.21) 

vlo vl5 



By considering additional Cartan generators beyond E and X, more embeddings of 50(10) can be found. This is discussed further in 
the appendix. 
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according to Eq. (|4.7p . The {X, E) basis is the more convenient for solving the Euler-Lagrange equations, because E 
and X are orthogonal as per Eq. (|4.6p . The {E,E') basis, however, is the simplest one for thinking about the two 
embeddings. 

The VEVs we want for the boundary conditions are 

{fE,.fE')^v{l,0), (4.22) 

which corresponds to Eq 50(10) U{1)e- The other VEV is 

(/b,/b') = -«(0,1) (4.23) 

which gives Eq 50(10)'® U{1)e'- The relative minus sign between Eqs. I|4.22p and (|4.23p comes from the breaking 
of a 

X^-X (4.24) 

discrete symmetry we shall impose on the Higgs potential, and it is crucial for two reasons. First, the spontaneous 
Z2 breaking will ensure that our domain walls are topologically non-trivial. Second, it leads to a remarkable outcome 
for fermion zero-mode localisation, to be explained in the next subsection. 
In terms of the {X, E) basis, these same VEVs are 

{fxjE)^v{0,l) and v(^^,-^, (4.25) 

respectively. 

We now need to find a Higgs potential with these two VEVs as degenerate global minima. The Higgs potential is 
constructed out of adjoint invariants, which according to Eqs. (|4.19p and (|4.18p are 

/„ = Tr(x") = Tr(V"il"^ • • -l^'O/aJa. •••/«„• (4.26) 



They are simply the nth order Casimir invariants. According to Refs. 38|, |39|, the independent invariants are 



/2, h, h, h, I9, /i2, (4.27) 

which immediately has an interesting consequence: the fact that I^^g are nonzero means that the discrete Z2 of 
Ea. (|4.24p is not a subgroup of Eq, because imposing it eliminates the otherwise present odd-power invariants. 
It is sensible to truncate the Higgs potential at order-six: 

V = -\ih + Mhf - 2304K/6 + ^A3(/2)^ (4.28) 

where some peculiar numbers and signs have been inserted for later convenience. In the extra-dimensional setting, 
field-theoretic models must generally be considered as effective theories valid below an ultraviolet cutoff scale A, 
because they are almost inevitably non-renormalisable. In writing down a Higgs potential, one simply adds terms of 
ever higher mass-dimension and truncates appropriately, given that the higher the mass-dimension the more suppressed 
it should become. For the Eq application, it is not helpful to truncate at fourth order, because the only fourth-order 
invariant is and I2 is invariant under an accidental 0(78) symmetry. The presence of Iq reduces the symmetry 
of the Higgs potential to Eq (presumably), and eliminates a pseudo-Goldstonc boson issue. ^ 

Equation (|4.28p is a complicated sextic in seventy-eight fields. But to perform the global minimisation analysis, 
one can always transform any VEV pattern to a standard form given by linear combinations of just the six generators 



^ An alternative is to truncate the classical theory at fourth order, but to add a Coleman- Weinberg potential generated through quantum 
corrections that explicitly break the 0{78) |39ll . 
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in the Cartan subalgebra of Eq. This produces a still quite complicated sextic in six fields. To make our discussion 
as simple as possible, in the main body of the paper we shall further truncate to just the two Cartan subalgebra 
generators of interest, and use Eq. (|4.20p . We extend the global minimisation analysis to all six fields in the appendix. 
With just fE,x 7^ 0, the nth-order invariant simplifies to 

^«^E( " J Tr(i?"-'=X'=)/r'/i. (4.29) 
k=o \ / 

The traces can be worked out by hand, because we know the matrix representations of E and X from the branching 
rules in Eq. (|4.4p . We obtain 

h - \{fl + fx), (4.30) 

/6 = ^ (/I + 5/l/i + 7/l/i - ^jEfx + • (4.31) 
To understand the extrema of Eq. (|4.28p , it is helpful to use the polar decomposition 

fE^r cos 9, fx =r sin 9. (4.32) 

The VEVs of Eq. are then 

(10,+) : 9 = 0, (4.33) 



(10',-): cos6l = -i, sin6'= (4.34) 

The notation (10, +) means that the VEV of Eq. (|4.33p induces Eq — > 50(10) U{1)e, and we have (arbitrarily) 
assigned it a positive Z2 "parity" which signals that it lies in G/H rather than {G/H)^. Similarly, (10',—) means 
Eq — > 50(10)' O U{1)e' and it lies in {G/H)z. There is another pair, with the opposite Z2 parities: 

(10,-) : 9 = TT, (4.35) 



(10',+) : cos6l = i, sin6' = (4.36) 

The topological CoS domain wall connects (10,+) and (10',—), accompanied by a CoS anti-domain-wall connecting 
(10,—) and (10',+). The topological non-CoS domain walls connect (10,+) with (10,—) [breaking Eq to 50(10) (8) 
U{1)e at aU y], and (10', +) with (10', — ) [breaking Eq to 50(10)' (g) U{1)e' at aU y]. There are also nontopological 
configurations: (i) (10,+) connected to (10',+), and (ii) (10,—) connected to (10',+), which are both CoS-like. 

Figure [2] display the invariants — /g/r^, — /g/r^, — (/s)^/?'^*' and — /i2/?'^^ as functions of 9. They show a remarkably 
similar structure. It is evident that the global minima for all four invariants are precisely the four VEVs of Eqs. (j4.33^ 
I4.36P . It is clear from this that choosing to truncate at the sextic level, as in Eq. (|4.28p , does not sacrifice much in terms 
of generality. We can be confident that our simplified potential leads to solutions whose qualitative characteristics 
would be retained were a wider class of higher-order potentials considered. In addition, the appendix shows that there 
are no deeper minima than (10, ±) and (10', ±) in the whole six-dimensional Cartan domain. 

The sign in front of Iq must be negative to achieve the desired extrema as minima rather than maxima. The other 
terms in the Higgs potential, Eq. (j4.28p . are independent of 9, depending only on the radial function r. Hence, it is 
clear that Eq. (|4.28p has the global minima we require. Figure [3] shows a contour plot of the Higgs potential for a 
certain parameter choice illustrating this conclusion. It is important to realise that although the minima (10, +) and 
(10', +) [similarly (10,—) and (10',—)] look as though they are disconnected by Eq, this is just an illusion created 
by only plotting the two-dimensional {fx, fs) slice through the 78-dimensional adjoint representation space. Minima 
with opposite parities are definitely disconnected from each other. 



While it is certainly true that the Z2 is not a subgroup of Eg, so that in general x E'nd ~X ^''c disconnected from each other, one may 
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FIG. 2: The invariants — Je/r®, —Ig/r^, —{Ir,)'^/r^" and ~Ii2/r^^ as functions of 9. Note the remarkable similarity. The global 
minima for all four, reading from left to right, are: (10,—), (10', -|-), (10, -|-), (10',—) and then (10,—) again (see text for 
explanation of this notation) . 




-0.8 -0.4 0.4 0.8 

fE 

FIG. 3: Contour plot of the Higgs potential as a function of the two field components Je and fx- The parameters used are 
n = 0.8, Ai = 1.0, A2 = 0, A3 — 22.0. The darkest regions are the global minima in the order (10,+), (10',—), (10,—) and 
(10', -|-) reading anticlockwise from the rightmost minimum. The light area near the origin is a local maximum. 



worry that there is nevertheless an Ea transformation that takes the specific configuration (10,+) to (10,—). However, this is not the 
case. Explicit calculation of the i?g invariant I5 reveals a nonzero f'^ term. Hence, /e ^fs must be outside of Ee- It does not 
matter that I5 has been omitted from the Higgs potential, as it is a purely group-theoretic argument. 
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-20 -15 -10 -5 5 10 15 20 -20 -15 -10 -5 5 10 15 20 

extra-dimensional coordinate y extra-dimensional coordinate y 

FIG. 4: Clash-of-symmetries domain wall solutions interpolating between (10, -|-) aX y = —oo and (10',—) at j/ = -foo. The 
parameters used in the left plot are k = 0.2, Ai = 1.5, A2 = 0, A3 = 22.0; those in the right plot are k = 0.8, Ai = 1.0, A2 = 0, 
A3 = 22.0. 



In the examples presented below, we further simplify the Higgs potential by setting A2 = as this term does not 
play an important role. It is then easy to show that at the degenerate minima, 

so we must take A3 > 22k, and that the value of V at the minima is 
The latter must be subtracted from the potential 

V^V- K.i„ (4.39) 

to produce finite energy-densities for the domain wall configurations. Figure[3]is a contour plot of the potential energy 
showing the four degenerate global minima. 

Having understood the global minima, we may now solve the Euler-Lagrange equations 

dV dV 

fx{y) = 7T^, my) = jrr (^.40) 

OJX OJE 

using those VEVs as boundary conditions. Numerical solutions for CoS domain walls interpolating between (10, -|~) 
at y = — oo and (10', —) aX y — -\-oo with two different parameter choices are displayed in Figure 21 

Figure [5] depicts non-CoS domain wall solutions for the same parameter choices. The function fx is zero, while 
fE interpolates between v and — u in archetypal kink fashion. This means that the non-CoS configurations feel the 
large potential-energy maximum at fx = /e = 0, while the CoS configuration "skirts around" that central maximum. 
This immediately implies that the CoS solutions have lower energy density than the non-CoS solutions. Although 
they are in the same topological class, the CoS domain walls are stable while the non-CoS domain walls arc unstable. 
Figure [6] shows a three-dimensional plot of the potential and where the two DW configurations sit with respect to the 
topography. There is a tall maximum at the origin, and a corrugated valley encircling it, with four low points at the 
VEVs. Figure [7] compares the energy densities of CoS and non-CoS domain walls. 



C. Fermion zero- mode localisation 



The CoS Eq domain wall solutions described above are a good starting point for the creation of domain-wall- 
brane models featuring S'C/(5)-invariant effective 3 -I- 1-d theories for localised fields. To actually create such a 
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FIG. 5: Non-cl£ish-of-symmetries domain wall solutions interpolating between (10, +) a,t y = —oo and (10, —) at y = +oo. The 
parameters used in the left plot are k = 0.2, Ai = 1.5, A2 = 0, A3 = 22.0; those in the right plot are k = 0.8, Ai = 1.0, A2 = 0, 
A3 = 22.0. 




FIG. 6: Three-dimensional plot of the Higgs potential as a function of the two field components fs and fx- The white lines 
show the clash-of-symmctrics domain wall (topmost) and the non-CoS domain wall (bottommost). The parameters used are 
K = 0.8, Ai = 1.0, A2 = 0, A3 = 22.0. 



model, fermions, additional Higgs bosons and gravitons have to be added. In this subsection, we demonstrate that 
a phenomenologically-acceptable fermion localisation pattern is obtained using the simplest possible mechanism. We 
explain why this is a remarkable result. 

We simply Yukawa-couple a five-dimensional fermion multiplet in the 27 of Eq, 

* ~ 27, (4.41) 

to the adjoint scalar, as per 

Cy = -h^x^- (4.42) 

We now substitute in the background CoS DW configuration for x and solve the resulting Dirac equations, which 
take the form 



iT^dM^™{x'^,y) - h[fxiy)X + fE{y)E]^Sf^'''^\x>^,y) = 0. 



(4.43) 
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FIG. 7: The difference in energy densities between the non-CoS and CoS domain wall solutions, Snon-CoS — EcoS- We have 
numerically scanned through the parameter space with < k < 1 along the horizontal axis, and each successive curve represents 
a different Ai, beginning at Ai = 0.05 at the bottom and increasing in steps of 0.05 to Ai — 1.5 at the top. The energy difference 
is always positive, so the CoS domain wall has a lower energy. We set A2 = for simplicity. 
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FIG. 8: Clash-of-symmetries fermion localising profiles interpolating between (10, +) at y — —00 and (10', —) a,t y — +00. The 
parameters used in the left plot are k = 0.2, Ai = 1.5, A2 = 0, Aa = 22.0; those in the right plot are k — 0.8, Ai — 1.0, A2 — 0, 
A3 — 22.0. The top to bottom order of the SU{5) fermion multiplets in the box on the right matches the order in Eq. (|4.48|l . 



The notation \I>(^'^) signifies the component of the 27 with the specified {X,E) charges, as given in Eq. (14. 4p . The 
various components couple to different background field configurations, 

b^''^''^ (y) ^ fx{y)X + fE{y)E, (4.44) 

given by the DW configuration and the charges. 

The Dirac matrices are = (7*^, —475). We search for separated- variable solutions 

^-(2:^2/) =F(y)V'(x'^), (4.45) 

demanding that have definite chirality, jr^ip = ±"0, and obey the 3 + 1-d massless Dirac equation, i^^d^^p — 0. The 
solution for a profile is well known: 

^^(^■^)(y) =iV(^'^)e-"/''^'"'"'(^')'*''', (4.46) 

where is a normalisation factor. For the profile to represent localisation, it must be square-integrable with respect 
to y. For this to happen, b'^'^'^'^ must pass through zero. If so and it is an increasing function of y (kink-like), 
then a left- (right-) handed zero-mode occurs for h > 0{h < 0). If it passes through zero as a decreasing function 
(antikink-like), then a left- (right-) handed zero-mode occurs for h < 0{h > 0). 
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Figure [8] show the kink-like functions for the two parameter choices we have been using as examples. Let us 

take h to be negative: 

h<0. (4.47) 

Using the notation £)(4\/T5-A:,12£;) once again for the fermion multiplets, the following displays these functions and states 
the localisation outcome, which is either "localised as left-handed (LH) zero-mode" or "localised as right-handed (RH) 
zero-mode" or "delocalised": 

1^°'''^ : ^/iJ, localised LH 
5(2-2). ^f^J^f^_^f^y localised RH 

5*(-2'-2) : -i [-^fx + IfE^ delocalised 
l(-5.i): U-^fx + ^fEY localised LH 



5*(3,i): l(^^f^ + lf^y delocalised 

w\Zli] ■ \ (";^^^ + J^^) ' localised LH. (4.48) 

The two 5*'s are delocalised because the associated field never goes through zero. The 5 and the 10 are localised at 
y = Q with opposite chiralities because their background fields are kink-like and antikink-like, respectively. The two 
singlets are localised at nonzero y values, so the overall spectrum is "split". 

This is a remarkable outcome for two reasons. First, because the 5 is localised RH, it is equivalent to a LH-localised 
5*. Thus the localised spectrum consists of LH zero-modes in the SU{5) representation 

5*®10®1©1, (4.49) 

in other words one standard family plus two singlet neutrinos. Second, apart from the extra singlet, all the exotic 
fermions in the 27 of Eq are delocalised and thus do not feature in the effective 3 -I- 1-d theory on the brane. These 
benign outcomes depend crucially on the boundary condition choice embodied by the CoS domain wall solution, 
including the Z2 minus sign. 

Finally, there is an amusing aspect to this spectrum. It resembles a usual 5*0(10) family plus an extra singlet. 
However, the LH 5*, which is obtained from a 4 + 1-d 5, and the 10 do not come from a 16 of either S'O(IO) or 
50(10)'. 



V. CONCLUSION 

We find it extremely encouraging that, in the Eq context, the clash-of-symmetries idea leads to good outcomes for 
both gauge-boson localisation (assuming the Dvali-Shifman mechanism works) and fermion localisation; that is the 
main point of this paper. 

In summary, we have established a general connection between the clash-of-symmetries mechanism for simultaneous 
brane-creation and internal-symmetry breaking with the Dvali-Shifman mechanism for gauge boson localisation. The 
two together provide a strong basis upon which to construct realistic domain-wall-brane models. These models should 
be compatible with type-2 Randall-Sundrum graviton localisation (see Ref. 16] for a CoS-style toy model featuring a 
background warped metric). 

More specifically, we have found a domain wall solution in an Eq adjoint-Higgs model that produces an 5[/(5)(8)C/(l)2 
symmetry on the wall itself. In one half of the bulk, the symmetry is enhanced to 50(10) ^ U{1), while in the other 
half of the bulk the enhancement is to 50(10)'® ?7(1)'. The unprimed and primed groups are differently-embedded but 
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isomorphic subgroups of Eq. Because the brane-S'[/(5) is contained in both 50(10) and 50(10)', the Dvah-Shifman 
locahsation of its gauge bosons follows. The simplest possible mechanism for fermion zero- mode localisation produces 
a realistic spectrum, an outcome that depends on the generic features of our domain wall configuration. 

To complete a realistic model, one needs to add gravity (which is expected to be straightforward) and to arrange for 
the additional spontaneous symmetry breaking cascade SU{5)(»U{1)^ SU{3)(»SU{2)(^U{1)y SU{3)(S)U{1)q. To 
achieve the latter, suitable additional Higgs multiplets need to be introduced, and their background field configurations 
have to be nonzero inside the domain wall to trigger the additional spontaneous symmetry breaking. An example of 



this kind of dynamical structure is described in Ref. [13|, where the dominant background domain- wall configuration 
breaks SU{5) to SU{3) (g) SU{2) (g) U{1), and then an additional Higgs field induces electroweak symmetry breaking 
inside the wall. 
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APPENDIX A: FULL MINIMISATION ANALYSIS 

We do this in two steps. We first extend the Higgs potential minimisation analysis by adding a third adjoint 
component, associated with the Cartan sub- algebra generator identified as weak hypercharge Y. This is useful 
because the result can be graphically visualised, and it reveals a third embedding of 50(10) that is related to the 
two embeddings used in the main body of the text. In the second step, we report on a numerical study of the whole 
six-dimensional Cartan subspace. 

The truncated multiplct is thus first increased to 

fEE + fxX + fvY, (Al) 

where 

E = ^diag(4, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), 

X = ^^diag(0, 2, 2, 2, 2, 2, -2, -2, -2, -2, -2, -5, 3, 3, 3, 3, 3, -1, -1, -1, -1, -1, -1,-1,-1, -1,-1), 
4vl5 

V - -^H- fn 1 ^ ^ ^ ^ ^^n 222111111 4 4 4 \ 

" 27TU ^""^V ' 3'"3'"3'"^'"^'3'3' '"^'"^'3'3'3'3'3'3'3'3'3'"3 ^ ^ 

have been normalised as per 

Tt{E^) = Tr{X^) = Tr(y2) = i/2, Tt{EX) = TiiEY) = Tt{XY) = 0. (A3) 
The sextic invariant is 

h = TtHJeE + fxX + fYY)% (A4) 
To visualise its structure, we go to a spherical-polar decomposition 

/e = r sinipcosO, fx = r siiKpsinO, fy = r cos 4>, (A5) 
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which produces 



710 cos^ (/) + 30 cos^ (/) (^51 + 4 cos 26* - 4^15 sin 26*) sin^ 



518400 

60\/2 cos^ (/) sin 61 ( 2^3 + 3 Vs cos 261 - V5 sin 29 I sin 
45cos2 6 f-34 + 2 cos 26* + 7 cos 46* - 2%/l5sin26' + Vl5sin46') sin'' 



+ ^ (^440 + 15 cos 26* + 84 cos 40 + 11 cos 661- 15\/T5 sin26' + 12Vl5sin40 - 3VT5sin66l) sin^ </> 

+ 144cos^0(yiOcos6'sin(/)- ^/6sin6lsin(^!)^^ (A6) 

Figure ini plots — /g/r^ as a function oi 9 and 0. The {E,X) plane is the line 4> = 7''/2, along which the VEVs (10, ±) 
and (10', ±) can be seen. Degenerate with them are two more VEVs with fy 7^ 0, located at 

arccos ( 7=] (A7) 

\ ^ V V VToJ 

and 

' ^ ' (A8) 

The first VEV corresponds to a nonzero value for the adjoint component associated with the generator —E", where 






E" = -E - -J-X + ^=Y. (A9) 



As the notation suggests, this minimum breaks -Eg to yet a third differently-embedded subgroup which we can call 
50(10)" ® U{\)e", with negative Z2 parity: (10", -). The second VEV is just (10", +). 

The three groups 50(10), 50(10)' and 50(10)" share a common SU{2,)®SU{2) subgroup, but the 5/7(5) contained 
in 50(10) n 50(10)' is not a subgroup of 50(10)" (this is obvious, since E" contains an admixture of Y which is 
a generator of that SU(5)). One can imagine a domain- wall junction configuration that utilises all three of these 
embeddings for boundary conditions, but such a model would have a similar photon and Z-boson leakage problem as 
the warm-up example of Sec. IIIII 

The functions —Is/r^, —{h)'^/^'^^ and —li^lr^^ have exactly the same qualitative structure as —Iq/t^. Thus 
(10, ±), (10', ±) and (10", ±) will be the degenerate global minima for a large class of potentials. The quadratic 
invariant is 9, ^-independent, 

h = \{fl + fl + f^) (AlO) 

so one can simply add appropriate terms to the potential to ensure it is bounded from below, and to generate a 
definite value for r at the global minima. The positions of the global minima are determined entirely from the angular 
structure of the non-isotropic terms. 

One can extend the analysis to all six Cartan components using a six-dimensional hyperspherical polar decompo- 
sition. The six fields are represented by one modulus, r, four zenith angles < '/>i.2,3,4 < tt and one azimuthal angle 
— TT < < TT. Because the group theoretic character of an extremum is determined entirely from the angular structure 
of the invariants, a numerical study can readily be performed on the finite domain (01.2,3,4, 9). This study confirmed 
that the Eq 50(10) ® U{\) VEVs are the global minima of — /6,8,i2 and — (15)^. The additional field dimensions 
simply revealed more degenerate vacua, corresponding to extra embeddings of 50(10) in Eq. These new embeddings 
must correspond, physically speaking, to choosing different 5[/(3) ® SU{2) subgroups for colour and isospin. The 
total number of Eq 50(10) (E) U{1) extrema was found to be 54, consisting of 27 Z2-related pairs. This implies 
that, overall, there are 27 embeddings of 50(10) (E) U{1) in E^. Though we shall not display the results here, we 
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-Ji -nil %I2 % 

9 



FIG. 9; Contour plot of -Je/r^ as a function of and <f). The degenerate global minima breaking E(, to various 50(10) 8 U{1) 
subgroups are marked with + signs, while the local maxima are indictaed with x signs. The row of minima along = 7r/2 
correspond to the global minima displayed in Fig. |2] 

have analytical expressions for the 27 linear combinations of Cartan generators that correspond to these VEVs. In 
the breakdown 27 ^ 1 ® 10 ® 16, these 27 linear combinations turn out to be correlated with the choice of which 
component to assign as the 5*0(10) singlet in the decomposition. A deeper reason for the number 27 is perhaps the 
following: according to the 5?7(3) ® SU{3) (X) SU{3) maximal subgroup of Eq, there are three independent choices 
for the colour group. The weak-isospin group can then be selected as the /— , U— or spin subgroup of either of 
the two remaining SU{3)^s. This gives 3 x 6 = 18 choices for SU{3) (8) SU{2) embeddings. According to our previous 
analysis, each SU{3) SU{2) is contained in the intersection of three different S'O(10)'s, which suggests there should 
be 18 X 3 = 54 embeddings of 50(10). However, recognising that 50(10) contains an SU{2) (g) SU{2) subgroup, we 
see that the correct number of independent embeddings is actually 54/2 — 27. 
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